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Abstract 



Many models of new physics, including versions of supersymmetry (SUSY), predict 
production of events with low missing transverse energy, electroweak gauge bosons, 
and many energetic final-state particles. The stealth SUSY model yields this signa- 
ture while conserving R-parity by means of a new hidden sector in which SUSY is 
approximately conserved. The results of a general search for new physics, with no re- 
quirement on missing transverse energy, in events with two photons and four or more 
hadronic jets are reported. The study is based on a sample of proton-proton collisions 
at = 7 TeV corresponding to 4.96 fb^^ of integrated luminosity collected with the 
CMS detector in 2011. Based on good agreement between the data and the standard 
model expectation, the data are used to determine model-independent cross-section 
limits and a limit on the squark mass in the framework of stealth SUSY. With this first 
study of its kind, squark masses less than 1430 GeV are excluded at the 95% confi- 
dence level. 
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Models of supersymmetry [[TJ m (SUSY) that conserve R-parity and include a neutral, weakly 
interacting lightest supersymmetric particle (LSP) garner much interest because they simulta- 
neously solve the hierarchy problem, allow unification of the fundamental interactions, and 
provide a candidate for dark matter. Many searches for SUSY focus on the presence of large 
missing transverse energy (E^^®) carried away by the LSP. It has been pointed out [3-5 1 that 
this approach neglects well motivated SUSY models that predict low E^^^ without special tun- 
ing of masses - models characterized by R-parity violation |6|, gauge mediated SUSY break- 
ing \7\, compressed spectra [8, 9], and hidden valleys IITOl . As the parameter space available 
for high-£^^*® SUSY is reduced by recent results from the CERN Large Hadron Collider (LHC), 
low-E^'^^ alternatives become more important to study. 

Without the discriminating power of large E^'^®, new techniques are required to reduce stan- 
dard model (SM) backgrounds. In general, searches in final states with many energetic par- 
ticles, including photons or leptons, are sensitive to the wide range of new physics models 
producing heavy particles with cascade decays to hadronic jets and electroweak gauge bosons. 
This range includes theories of extra dimensions l|TT| . heavy-flavor compositeness [12], little 
Higgs [ 13^. 14], and even SUSY. In particular, such searches are sensitive to the R-parity con- 
serving stealth SUSY model [5. 15\, which, as formulated in this study, produces events with 
low E^'^®, two photons, and more than five jets. 

The simplest stealth SUSY models assume low scale SUSY breaking and introduce a new hid- 
den sector of particles at the weak scale in which only a small amount of SUSY breaking occurs 
through interactions with SM fields, resulting in a hidden sector that is approximately super- 
symmetric and nearly mass degenerate hidden sector superpartners. In this framework, the 
standard LSP takes on a new role as the lightest "visible sector" SUSY particle (LVSP) which 
decays without violating R-parity into a lighter hidden sector SUSY particle. In the subsequent 
decay of this particle to its SM partner and the true LSP, the near mass degeneracy leaves little 
phase space for the true LSP to carry momentum. In this way, stealth SUSY models naturally 
produce signatures of low E^'^® without special tuning of masses. 

In this Letter we describe a search for new phenomena (NP) in events with at least two photons 
and four or more jets produced in proton-proton collisions at ^/s = 7TeV. The analysis is 
based on a data sample corresponding to 4.96 ± 0.11 fb^^ of integrated luminosity [16] collected 
in 2011 with the Compact Muon Solenoid (CMS) detector at the LHC. We perform a model- 
independent search for an excess of events with respect to the SM expectation, and we interpret 
this search in the framework of stealth SUSY. To minimize the impact of SM backgrounds on the 
sensitivity to NP, event excesses are determined as functions of final-state jet multiplicity and 
the scalar sum of the transverse momentum pj of jets, photons, and E^^^ (for E^^^ > 20 GeV) 
in each event. This sum is referred to as Sj; the definitions of the objects included in the sum 
are described below. 

We consider a model that includes degenerate light squarks (q), a "bino-like" LVSP x^, a gluino 
(g) with mass of 1500 GeV, and a hidden sector containing a singlet state S and its fermionic 
"singlino" superpartner S. The model is similar to squark-antisquark (q q*) production with 
the decay q — > q^ described in Ref. [17] with two differences: the addition of the hidden sector 
and the participation of g in the production mechanism. After initial production of two ^ and 
jets, each ^ decays into the hidden sector producing a photon and S, which subsequently 
decays to S and a gravitino, S — > SG. The S state is even under R-parity and decays to jets via 
S— >gg. The resulting G LSP has small momentum because the hidden sector superpartners (S, 
S) are nearly mass degenerate, so the final state tends to have low E™*®. The decay of a squark 
in this model is shown in Fig. [l] 
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Figure 1: Decay of a squark in stealth SUSY. See text for a complete description. 

The model is characterized by the masses of the particles in the decay chain. We consider a 
range of squark masses Mq from 400 to 2000 GeV in steps of 100 GeV. We make the following 

assumptions inspired by benchmark points described in Ref. ||5l: S and S masses of 100 GeV 
and 90 GeV, respectively; with mass equal to ^Mg; and branching fractions of unity for the 
decays of x^, S, and S described above. The production cross section for this process (and its 
uncertainty) is calculated as a function of Mq at next-to-leading order (NLO) accuracy includ- 
ing the resummation of soft gluon emission at next-to-leading logarithmic (NLL) accuracy as 
described in Ref. llT8ti23ll . 

A detailed description of the CMS detector can be found in Ref. H24| . The CMS coordinate 
system is right-handed with the origin at the center of the detector, the x-axis directed toward 
the center of the LHC ring, and the y-axis directed upward; (p is the azimuthal angle, 9 is the 
polar angle, and tj = — ln[tan(0/2)] is the pseudorapidity. The central feature of the CMS 
apparatus is a superconducting solenoid of 6 m inner diameter that surrounds a silicon pixel 
and strip tracker, covering the region \rj\ < 2.5, as well as a lead-tungstate crystal electromag- 
netic calorimeter (ECAL) and a brass/ scintillator hadron calorimeter (HCAL), both covering 
\r]\ < 3. Muons are detected by gas-ionization detectors embedded in the steel flux return yoke 
covering the range |?/| < 2.4. 

The triggers, event reconstruction methods, and selection criteria for photons and jets are iden- 
tical to those used in the CMS search for SUSY in events with photons, jets, and E^'^® |j25j. 
Events are recorded with the CMS two-level trigger system requiring the presence of one pho- 
ton with transverse energy (£x) greater than 36 GeV and a second with Ej > 22 GeV. To 
suppress jets giving rise to photon candidates, these triggers require the latter to be isolated 
from other activity in the tracker, ECAL, and HCAL. As instantaneous luminosity increased 
throughout 2011, isolation requirements were gradually changed to keep the trigger rate ap- 
proximately constant, but the isolation in the trigger is always less restrictive than offline re- 
quirements described below. 

Photon candidates are reconstructed from clusters of energy in the ECAL barrel with | ?/ 1 < 1 .44. 
Candidate events are required to have a leading photon with Ej > 40 GeV and an additional 
photon with Ej > 25 GeV; at these thresholds the triggers are more than 99% efficient. We 
require the ECAL cluster shape to be consistent with that expected for photons, and the energy 
detected in HCAL in the direction of the photon shower to not exceed 5% of the ECAL energy. 
We ensure isolation from other activity in the event by requiring that the scalar Ej sum of 
tracks and calorimeter deposits within AR = ^ {Ar]Y + {A(pY = 0-3 of the photon candidate's 
direction be less than 6 GeV after correcting for contributions from the products of additional 
collisions in the event (pile-up) and the candidate itself. These criteria efficiently select both 
photons and electrons; candidates that cannot be matched to hit patterns in the pixel detector 
are considered photons. 
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Jets are reconstructed with the particle-flow algorithm l|26l , which simultaneously reconstructs 
all particles produced in a collision based on information from all detector subsystems and 
identifies each as a charged or neutral hadron, photon, muon, or electron. All of these particles 
are clustered into jets with the anti-fcj clustering algorithm [27J with radius parameter of 0.5. 
To remove jets arising from potential instrumental and non-collision backgrounds, we require 
the fraction of jet energy coming from charged and neutral electromagnetic deposits to be less 
than 0.99, the neutral hadron fraction to be less than 0.99, and the charged hadron fraction to 
be greater than zero. The jet energy and momentum are corrected for the nonlinear response 
of the calorimeter and the effects of pile-up. Jets are required to have pj > 20 GeV, < 2.4, 
and to be isolated from photon candidates by AR > 0.5. The fraction of selected events that 
includes an isolated lepton in addition to photons and jets is less than 0.1%. 

We determine the fraction of stealth SUSY events that pass these reconstruction procedures 
and selection criteria (the "acceptance") and the leading order plus leading logarithm (LO) 
cross section for stealth SUSY using the PYTHIA 6.424 event generator (281, the D6T underlying 
event tune II29II , the CTEQ6L1 [30 1 parton distribution functions (PDFs), and a full simulation 
of the CMS detector based on GEANT4 [31 1. The acceptance when including a requirement 
for five or more jets rises monotonically from 17% to 35% for Mq from 400 to 1100 GeV and falls 
monotonically from 35% to 25% for Mq from 1100 to 2000 GeV as additional hadronic activity 
makes it less likely that photons will be isolated. The acceptance when requiring exactly four 
jets is 12-13% (1-2%) of the >5-jet acceptance for Mq < 1500 GeV (Mq > 1500 GeV). The LO 
acceptance for each of the dominant subprocesses (pp — > q q/ q q*/ q g/ and g g) is weighted with 
the appropriate NLO/LO K-factor; the corrected acceptance is within 5% of the LO acceptance 
at all Mq. 

Events satisfying the selection criteria come mainly from direct production of two photons in 
the initial parton-parton hard scatter with radiation via quantum chromodynamics producing 
the additional jets. These events are divided into three samples: a signal-rich "search sample" 
comprising events with four or more jets and Sj > 700 GeV, a signal-depleted "jet multiplic- 
ity sideband" (JMSB) composed of events with two or three jets and Sj > 600 GeV, and a 
signal-depleted "Sj sideband" composed of events with four or more jets and 600 GeV < Sj < 
700 GeV. We examine the search sample for evidence of NP in the form of an excess of events 
over the SM expectation. The SM expectation is estimated from the data based on the observa- 
tion that the shape of the Sj spectrum of the SM background is independent of jet multiplicity: 
we take the Sj shape from the JMSB and the normalization from the Sj sideband. 

This method of background estimation was first used in the CMS search for black holes p2ll33| , 
in which the jet-multiplicity invariance of the Sj shape in jet dominated events was demon- 
strated. We confirm that this invariance holds for events with photons in addition to jets using 
a data sample of events with one photon and two or more jets (photon+jets). The event selec- 
tion criteria for this sample are the same as those described above except for changes required 
by differences in the trigger and event topology: we require that events include a single photon 
with Ex > 80 GeV and Hj > 450 GeV, where Hj is defined as the scalar sum of the pj of jets 
with Pj > 40 GeV and \r]\ < 3.0. 

In Fig. |2]we compare the Sj spectra for five subsamples of this photon+jets dataset character- 
ized by jet multiplicity: the 2-jet, 3-jet, 4-jet, >5-jet, and >3-jet samples. We show these spectra 
(area-normalized for Sj > 800 GeV) along with the n-jet/2-jet ratio where n = 3, 4, >5, and >3. 
The shape of each ratio is consistent with a flat line within the statistical uncertainty. As an ex- 
ample, the fit of the >3-jet/2-jet ratio from Fig.|2j3 with the function a + hx where x = Sj/7Te\f 
yields a = 1.0 ± 0.2 and h = —0.3 ± 1.1 (statistical uncertainty). A slope of this size would 
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affect the expected background rate by 6%, which is negligible compared to the systematic 
uncertainty from other sources (described below). 




Figure 2: Sj spectra from the photon+jets data control sample, area-normalized for Sj > 
800 GeV. We show (a) the spectra for events with two, three, four, and five or more jets along 
with the n-jet/2-jet ratio where n = 3, 4, or >5; and (b) the spectra for events with two and three 
or more jets along with the >3-jet/2-jet ratio. 



As introduced above, we model the shape of the Sj distribution for the background by fitting 
the JMSB data with the function l/x"^, where x is defined above and a is a free parameter. We 
find a = 5.02 ± 0.32 from the fit of the entire JMSB; in the 2-jet (3-jet) subsample of the JMSB 
we find a best fit value of a = 5.00 ± 0.45 (5.03 ± 0.45). This background shape is normalized 
for use in the search sample with data in the Sj sideband. 

We compare this background prediction to the observed Sj spectrum in the search sample in 
Fig. |3] along with the systematic uncertainty on the prediction (described below). The data are 
in good agreement with the background expectation. The probability for the most significant 
excess (at 700 GeV < Sj < 750 GeV) to have local significance as high or higher than that 
observed is 0.08. 

Supported by the good agreement of the data with the background expectation, we use the data 
to compute Nf^{Sf^), the upper limit on the number of events exceeding the SM expectation 
as a function of lower Sj threshold. Limits are computed at the 95% confidence level (CL) with 
the modified frequentist CLg method [34, 35J based on a profile likelihood ratio test statistic 
constructed from the Poisson probability for the observed number of events given the expecta- 
tions for background and signal. For the model-independent cross-section limits, we separately 
compute Ng™(SJp™) for two categories of jet multiplicity, >4-jet and >5-jet. For stealth SUSY 
limits, we compute Ng^(S^"'^) using exclusive 4-jet and >5-jet jet-multiplicity bins combined 
as described in Ref. II36I : the inclusion of the 4-jet bin improves the expected cross-section limit 
by 7.5% at Mq=1400 GeV. 

Using the measured integrated luminosity, we convert Ng^(S^"^) into a model-independent 
limit on the product of the acceptance and cross section of any new process that may be present 
in addition to the SM (the "effective NP cross section"). Using the integrated luminosity and 
the signal acceptance, we convert Ng™(S^'") into a limit on the stealth SUSY cross section as a 
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Figure 3: Observed Sj spectrum, background expectation with systematic uncertainty, and 
predicted signal for a squark mass of 900 GeV in events with four or more jets. 

function of Mq. For each value of Mq, we use the value of that maximizes the sensitivity to 
the stealth SUSY signal accounting for systematic uncertainty on the background expectation; 
these Sf^"" values range from 800 to 2400 GeV for Mq from 400 to 2000 GeV. 

Systematic uncertainties on the expected numbers of background events, the signal acceptance, 
and the integrated luminosity are treated as nuisance parameters with a Gamma function prior 
distribution for the uncertainty related to the normalization of the background prediction and 
log-normal prior distributions for the other uncertainties. The theoretical uncertainty on the 
predicted cross section does not enter the cross-section limit, but is used in determining the 
limit on Mq. 

We estimate that the following sources of uncertainty affect our knowledge of the expected 
background at the stated level: statistical uncertainty from the background normalization method 
(15%), statistical uncertainty on the background shape (10-55%), and the choice of background 
function (5-100%) (ranges indicate dependence on Mq and SJJ^"^). We use the change in back- 
ground expectation obtained by varying the shape parameter by ±1 standard deviation (SD) 
as the statistical uncertainty on the background shape. We estimate the uncertainty related 
to the choice of background function by constructing background models from the alternate 
fit functions e"^ and l/x'^o+^i instead of the nominal l/x'^; we take the largest alternate- 
nominal difference in the background expectation as the uncertainty. The relative effect of these 
background uncertainties is large at high Sj where very few background events are expected. 

The following sources of uncertainty affect our knowledge of the expected numbers of signal 
events at the stated level: jet energy scale (3%), statistical uncertainty on signal acceptance 
from finite simulated samples (2%), and the measurement of integrated luminosity (2.2%) [16j. 
The effects on the signal acceptance of variations in pile-up and PDFs are less than 1%. The 
theoretical uncertainty on the predicted cross section related to PDFs, renormalization and fac- 
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torization scales, and ag variations is estimated to be 9-57% for Mq from 400 to 2000 GeV; the 
dominant source is the PDF uncertainty. 

Figure |4] provides the model-independent limit on the effective NP cross section as a function 
of for events with four or more jets and five or more jets. We show the observed limits, 
median expected limits, and a band corresponding to ±1 SD on the median expected limits. 
In Fig. |5] we show the same curves for the stealth SUSY cross-section limit as a function of 
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Figure 4: Model-independent limit at the 95% CL on the product of the acceptance and cross 
section for events with (a) four or more jets and (b) five or more jets. We show the observed 
limits, median expected limits, and a band corresponding to ±1 SD on the median expected 
limits. 

Mq. We also show the predicted NLO+NLL cross section with a band denoting ±1 SD from 
theoretical uncertainty. Based on the intersection of the cross-section limit and the — 1 SD edge 
of the predicted cross-section band, we exclude squarks with Mq<1430 
GeV for the stealth SUSY model described above. 

In summary, we perform a search for NP in events with two photons and four or more jets. The 
selection requirements are general and provide sensitivity to a broad range of NP phenomena. 
We observe no excess over the SM expectation in a data sample corresponding to 4.96 fb^^ 
of integrated luminosity collected in 2011. We determine model-independent limits on the 
effective NP cross section of 0.6 to 4fb depending on SJ^'". We also compute the limit on the 
stealth SUSY cross section as a function of Mq. Comparing this limit to the predicted stealth 
SUSY cross section, we exclude the production of squarks with Mq <1430 GeV. These excluded 
stealth SUSY points are characterized by low E^^^^ predicting 1.5 ± 0.1 (21.3 ± 1.6) events with 
E^'"" > 100 GeV for Mq of 1400 (800) GeV. This is the first limit on the parameters of the stealth 
SUSY model. 
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Figure 5: Stealth SUSY cross-section limit at the 95% CL as a function of Mq. We show the 
observed limit, median expected limit, and a band corresponding to ±1 SD on the median 
expected limits. We also show the predicted NLO+NLL cross section from stealth SUSY with a 
band denoting ±1 SD from theoretical uncertainty. 
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